Released momentum distribution of a Fermi gas in the BCS-BEC crossover 
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We develop a time-dependent mean- field theory to investigate the released momentum distribution 
and the released energy of an ultracold Fermi gas in the BCS-BEC crossover after the scattering 
length has been set to zero by a fast magnetic-field ramp. For a homogeneous gas we analyze the 
non-equilibrium dynamics of the system as a function of the interaction strength and of the ramp 
speed. For a trapped gas the theoretical predictions are compared with experimental results. 
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The study of the momentum distribution of an atomic 
gas in the quantum degenerate regime carries a wealth 
of information on the role played by interactions and on 
the existence of a superfluid order parameter. As an ex- 
ample, in a homogeneous Bose gas at T = 0, the mo- 
mentum distribution exhibits a singular behavior at 
small wave vectors ~ mc/2hk, which is determined 
by the excitation of phonons propagating with the speed 
of sound c and is a signature of Bose-Einstein conden- 
sation 1]. In a corresponding system of fermions with 
attractive interactions, the broadening of the Fermi sur- 
face is instead a consequence of the formation of pairs 
and of the presence of a superfluid gap 0- This latter 
effect becomes dramatic in the BCS-BEC crossover re- 
gion where the pairing gap is of the order of the Fermi 
energy of the system {3j . The second moment of the mo- 
mentum distribution defines the kinetic energy of the sys- 
tem: Ekin = J2k n\Ji 2 k 2 / (2m), where m is the mass of 
the atoms. This quantity, which also plays a central role 
in the many-body description of ultracold gases, is very 
sensitive to the large-fc behavior of n^. For interacting 
systems the dominant contribution to Ekin comes from 
short-range correlations, where the details of the inter- 
atomic potential are relevant. In the case of a zero-range 
potential it is well known that the momentum distribu- 
tion decreases like 1 /k for large momenta and the kinetic 
energy diverges in dimensionality greater than one. This 
unphysical divergence can be understood recalling that 
the zero-range approximation is only correct to describe 
the region of momenta k <C lM)j where r denotes the 
physical range of interactions This behavior of the 
kinetic energy is a general feature of quantum-degenerate 
gases, where interactions are well described by the s-wave 
scattering length a, holding both for fermions and bosons 
and for repulsive and attractive interactions p|. 

The physics of ultracold gases is characterized by a 
clear separation of energy scales: the energy scale asso- 
ciated with the two-body physics as fixed for example 
by h 2 /mrQ ~ lOmK, being ro ~ lOOao the typical in- 
teraction length of the Van der Waals potential, and the 
energy scale associated with the many-body physics as 
determined by the typical Fermi energy ep ~ \[iK . This 
separation of energy scales provides a very large range of 



timescales for which the dynamical process can be safely 
considered fast (diabatic) as the many-body dynamics is 
concerned and slow (adiabatic) with respect to the two- 
body dynamics. This feature is exploited in recent exper- 
iments aiming to measure the momentum distribution, 
that are based on the ballistic expansion of the cloud af- 
ter the scattering length has been quickly set to zero by 
a fast magnetic-field ramp 0, Q • These experiments give 
access to the released momentum distribution, which is a 
non-equilibrium quantity defined as the momentum dis- 
tribution of the system after the scattering length has 
been rapidly ramped to a = 0. Provided the timescale 
of the ramp satisfies the conditions given above, the re- 
leased momentum distribution does not depend on the 
detailed structure of the interatomic potential, being in 
this sense universal, but it does depend on the timescale 
of the ramping process. An important theoretical issue is 
to investigate the properties of the released momentum 
distribution and to understand what useful information 
about the system can be extracted from it. 

In this Letter we calculate the released momentum 
distribution of a Fermi gas at T = in the BCS-BEC 
crossover. We calculate the released momentum distribu- 
tion and its second moment for a homogenoeous system 
as a function of the interaction strength l/(kpa), where 
kp is the Fermi wave vector. For harmonically trapped 
systems we give an explicit prediction of the column in- 
tegrated released momentum distribution and of the re- 
leased energy for values of the interaction strength rang- 
ing from the BCS to the BEC regime and we compare 
our results with recently obtained experimental data 

We consider an unpolarized two-component Fermi gas 
with equal populations of the f and [ components: N-\ = 
N± = N/2, where N is the total number of particles. We 
determine the dynamical evolution of such a system start- 
ing from the equations of motion for the non-equilibrium 
density matrices of the j and j components interacting 
through the Hamiltonian 

H = J2j dx ^(x) ("^f ) ^(x) (1) 
+ / dxdxV|(x)V4(x')V(x,x')^ i (x , )V' T (x) 
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where a =f, J, labels the spins and V(x, x') is the inter- 
action potential to be specified later. The correlations 
are treated within a mean- field approach Q, where they 
are expressed in terms of the normal Gjvo-(x, x', t) — 
( , 0j.(x',t)'0 o -(x, t)) and the anomalous G^(x, x',i) = 
(ipl(yi', t)-0t( x 7 1)) density matrices. By neglecting the 
Hartree terms one obtains the following coupled equa- 
tions of motion 0] (from now on Gjvj (x, x', t) — 
G NL (x,x!,t) =Gjv'(x,x',i)) 



.. dGjv(x,x',i) 

zn, 

dt 
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Gjv(x,x',t) 



+ /" dx"[V r (x,x")-^(x",x')]G^(x,x",f)G^(x",x',t)(2) 
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dG^(x,x',i) 
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GA(x,x',t) 
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+7(x,x')G A (x,x',t) 

— / dx"V(x, x")Gjv(x', x", t)Gyi(x, x", t) 

-J dx"T/(x,x")Gjv(x,x",t)GA( x ",x',t) . (3) 

The short-range nature of the interaction potential 
V(x,x') can be properly described through the regular- 
ized pseudopotential V(r) = (inah 2 /m)S(r)(d / dr)r 0, 
with a the s-wave scattering length and r = |x — x'|. 
During the magnetic-field ramp the value of the scatter- 
ing length changes in time according to the relation 



(4) 



valid close to the Feshbach resonance. In the above ex- 
pression, cibg denotes the background scattering length, 
Bq and T the position and width of the resonance respec- 
tively, and B(t) is the instantenous value of the magnetic 
field. Under the dynamical conditions that we are con- 
sidering, where the non-equilibrium processes take place 
over a timescale adiabatic with respect to the two-particle 
problem and diabatic with respect to the many-particle 
system, the time evolution does not depend on the details 
of the short-range potential and the effect of interactions 
results in a boundary condition at short length scales 



(rG A (r,t))' 
rG A (r, t) 



1 



(•5) 



where the prime indicates the derivative with respect to 
r, which must be fulfilled at any time t. For small values 
of r, many-body effects in Eq. can be neglected and 
the boundary condition (JSJ) corresponds to the one of the 
two-body problem with the pseudopotential V(r), where 
G^fr, t) plays the role of the wave function for the relative 
motion. A similar argument is known to hold also in 
the electron gas, where one has the boundary condition 



Iff '( r ) 1 1 9{ r )\r=o — 1/clb for the pair correlation function 
g(r) (as is the Bohr radius), referred to as the Kimball 
relation Q. 

In the case of a homogeneous system and by us- 
ing the pseudopotential approximation for the inter- 
atomic potential V(x, x") and the boundary condition 
(0, Eqs. HJ-0 can be greatly simplified. For ex- 
ample, the term / dx"V{x,x")G N {x' ,x" ,t)G A (x,x" ,t) 
in Eq. © becomes (4nh 2 a/m,)G N (r, t)[(rG A )'] r = = 
— (Anh 2 /m)GN(r,t)[(rGA)]r=o, where we have also used 
the fact that the short ranged V(x, x") picks up the 
x" = x value of Gjv(x',x",t). After dealing with the 
other terms in a similar way, we obtain the simpler cou- 
pled equations for Gjv (r,t) = rGjv(r, i) and G A (r, t) = 
rG A {r,t) 



ih 



dG N (r,t) 
dt 



.^dG A {r,t) 

in 

dt 



8nh 2 
rn 

h 2 d 2 



3 (G A (r,t)[G A (t)] r=0 
G A (r, t) 



mdr 2 

—G N (r,t)[G A (t)]r=o 



(6) 



(7) 



with the boundary condition [(G a)' / G A ] r =o — — l/ffl(i). 
Notice that interaction effects only enter Eqs. ©-0 
through the boundary condition JSJ. We determine the 
initial conditions Gjv(r, i = 0) and G A (r,t = 0) of 
Eqs.|((jJ)-(|7J) from the mean-field gap and number equa- 
tions corresponding to the equilibrium state of the gas 
with the initial value of the scattering length a(0) 



(8) 
(9) 



Aith 2 a{0) 




e k ^{e k - Jl) 2 + A 2 



where = h 2 k 2 /2m, (x is the chemical potential, A the 
superfluid gap and n = n-f + the total particle density. 
The functions Gn and G A are then calculated from the 
Bogoliubov quasiparticle amplitudes 



2 -i 2 

u k = 1 - v k 



1 



V(e fc - M) 2 + A 2 



(10) 



as 



and 



G N (r,t = 0) 



2tt 2 



dkksin{kr)vl (11) 



G A (r,t = 0) = I dkksm(kr)u k v k 



(12) 



We solve the dynamic equations ®-0 with the initial 
conditions Qll [l -l|12 [l and a(t) given by Q from the initial 
time t = to the final time t = tj, where a(tf) = 0. The 
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FIG. 1: (color online) Released momentum distribution (solid 
lines) of a homogeneous gas at unitarity, l/(fci?a(0)) = 0, for 
a ramp rate of 2/^s/G. The large-fc behavior of «k weighted 
by fc 4 is shown in the inset, where the dotted line corresponds 
to the equilibrium asymptotic value (A/2ef) 2 . The initial 
equilibrium distribution is also shown (dashed lines). 




FIG. 2: (color online) Released energy as a function of the 
interaction strength for two values of the ramp rate. The en- 
ergy is normalized to the kinetic energy of the non-interacting 
gas E° km = 3e F /S. 



released momentum distribution is then calculated from 
the Fourier transform of Gn at the time t = tt 

n k (t = t f ) = J dre l * r G N (r,t = tf) . (13) 

The results for the homogeneous gas are shown in 
Figs. Ill'-'l In Fig. ^ we compare the equilibrium mo- 
mentum distribution n-^it = 0) in the unitary limit, 
l/(fc_pa(0)) = 0, with the corresponding released mo- 
mentum distribution (|13f) calculated for a magnetic-field 
ramp rate of 2/is/G. For values of A; < Uf the shape of 
the distribution does not change appreciably. The large- 
k tail is instead greatly suppressed (as is shown in the 



inset) and the second moment of the released n k is a 
convergent integral. Notice that the fast decaying tail 
of the released rik affects the normalization constant . In 
Fig. [21 we show the results of the released energy as a 
function of the initial interaction strength 1/ (fc^-a(O)) for 
two different values of the magnetic-field ramp rate. We 
notice that on the BCS side of the crossover, kpa(Q) < 0, 
the dependence on the ramp rate is weak, while on the 
BEC side, fcpa(O) > 0, a faster ramp produces a signif- 
icantly larger energy. In the BEC regime the system is 
in fact more sensitive to changes of the high-energy tail 
of the momentum distribution. Deep in the BCS regime, 
— l/(fci?a(0)) 3> 1, the released energy reduces to the ki- 
netic energy of the non- interacting gas E^. in — 3ef/8. 
In the opposite BEC regime, — l/(fcf a(0)) <C —1, many- 
body effects become less relevant and the released energy 
coincides with the one obtained from the dissociation of 
the molecular state (see Fig. 0) . 

In order to make quantitative comparison with the ex- 
periment, we now consider harmonically trapped systems 
confined by the external potential V ex t(r) = m(u} 2 x 2 + 
oj 2 y 2 + uj 2 z 2 )/2. Within the local density approxima- 
tion (LDA) we introduce the rescaled spatial variables 
x = xy/(jj x /iij, y = y^fujy/bj and z = x^f uj z /lu, so that 
the confining potential becomes isotropic in the new co- 
ordinates V ex t{v) = muj 2 R 2 /2, where ui = (uj x ujyUj z ) 1 /' i 
is the geometric average of the harmonic oscillator fre- 
quencies. For each spatial slice R = (x + x')/2, Eqs. (JSJ- 
© are solved for the local chemical potential ^ oca /(R) 
and the local density n(R), subject to the normaliza- 
tion Jd 3 Rn(R) = N and the local equilibrium con- 
dition fi = /i; oca /(R) + 14xt(R). Each slice is then 
evolved according to Eqs. ©-(0 with initial conditions 
Gjv(f, R, t = 0) and GA{r,R,t = 0), where r = x — x' 
is the relative coordinate. The released momentum dis- 
tribution is obtained from G;si{r,R,t) at the final time 
t — tf through the integral over the rescaled coordinates 
R and f , 

n{k,t = t f )= J d 3 R J d 3 re t]if G N (r,R,t = t f ) . (14) 

In Fig. |21 we compare the column integrated re- 
leased momentum distribution n(^J k 2 + k 2 ,tf) — 

J^° 00 'ik z n(k,tf), calculated from Eq. (|14fl . with the ex- 
perimental results obtained in Ref. [7| . The values of the 
interaction strength are l/(k F a(0)) =-0.66, and 0.59 
as in the experiments and the magnetic-field ramp rate 
is 2/is/G. The agreement is quite good for the momen- 
tum distribution on the BEC side of the resonance and 
in the unitary limit. On the BCS side of the resonance 
the experimental n(k) is more broadened because the 
Hartree mean-field term, which enhances the shrinking 
of the cloud due to attraction, is neglected in the calcula- 
tion. In Fig. 21 we show the results for the released energy 
of the inhomogeneous gas as a function of the interaction 
strength l/(k° F a(0)), where k F = (247V) 1 / 6 ^fmuj/h is the 
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FIG. 3: (color online) Column integrated released momen- 
tum distribution of a harmonically trapped gas. From top to 
bottom, the lines correspond to l/(k F a(0)) = —0.66 (green), 
l/(k° F a(0)) = (red) and l/(k° F a(0)) = 0.59 (black). The 
magnetic-field ramp rate is 2/is/G. The symbols correspond 
to the experimental results of Ref. 0. 
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FIG. 4: (color online) Released energy of a harmoni- 
cally trapped gas as a function of the interaction strength 
l/(k F a(0)) for a ramp rate of 2/xs/G (upper blu line). The 
lower (green) line is the corresponding result solving the two- 
body problem associated with the molecular state. The sym- 
bols are the experimental results from Ref. 0. The energy 
is normalized to the kinetic energy of the non-interacting gas 
E°in = 3ef/8- 



Fermi wave vector in the center of the trap correspond- 
ing to a non-interacting gas. Experimental results from 
Ref. Q and theoretical results obtained by solving the 
time-dependent Schrodinger equation for the molecular 



state (see 0) are also shown in Fig. The full many- 
body calculation reduces to the molecular two-body re- 
sult only in the deep BEC regime, — l/k F a(0) <C — 1, and 
agrees better with the experimental results. Considering 
that there are no adjustable parameters in the compari- 
son between theory and experiment, the agreement is re- 
markable over the whole crossover region. In the unitary 
limit the mean-field approach employed in the present 
study is known to overestimate the equilibrium energy 
per particle compared to more advanced quantum Monte 
Carlo calculations 10]. This might be the reason for 
the larger energy obtained around resonance compared 
to the observed one. Furthermore, on the BCS side of the 
resonance, the present approach neglects the mean-field 
Hartree term and we expect a faster convergence to the 
kinetic energy of the non interacting gas E® in = 3e^/8, 
where e° F = (hk F ) 2 /2m, as — l/fc^a(0) becomes large. A 
more sophisticated time-dependent theory is needed to 
improve the quantitative agreement with experiments. 

In conclusion, we have developed a time-dependent 
mean-field scheme which allows one to calculate the re- 
leased momentum distribution and the released energy 
of a Fermi gas in the BCS-BEC crossover if the scatter- 
ing length is set to zero using a fast magnetic-field ramp. 
Our analysis clarified that the dynamical effects due to 
the magnetic-field ramp suppress the high-energy tail of 
the momentum distribution. For harmonically trapped 
systems we compared our theoretical predictions with 
the recently obtained experimental results of Ref. 0. 
Qualitatively we reproduced the data well, but signifi- 
cant quantitative discrepancies are found. The underes- 
timate of the release energy of the theory compared to 
experiment on the deep BCS side is anticipated due to 
the neglection of the Hartree energies. On the deep BEC 
side, the underestimate may be due to the effect of finite 
temperature on the experimental data in that regime. 
The significant overestimate of the release energy at res- 
onance is more difficult to explain and may indicate the 
inadequacy of the Leggett ground state form for describ- 
ing the gas in the unitary limit. This would illustrate 
a need to include correlations beyond the pair correla- 
tions given in the Cooper pair picture in order to explain 
the observed experimental expansion energetics on reso- 
nance. 
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